Properties of holons in the Quantum Dimer Model 
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I introduce a doped two-dimensional quantum dimer model describing a doped Mott insulator 
and retaining the original Fermi statistics of the electrons. This model shows a rich phase diagram 
including a d-wave hole-pair unconventional superconductor at small enough doping and a bosonic 
superfluid at large doping. The hole kinetic energy is shown to favor binding of topological defects 
to the bare fermionic holons turning them into bosons, in agreement with arguments based on RVB 
wave-functions. Results are discussed in the context of cuprates superconductors. 
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The discovery of high temperature superconductivity 
in copper oxydes triggered unprecedented efforts to un- 
derstand the role of strong correlations in electronic sys- 
tems. In Anderson's milestone paper the "mother" 
correlated insulator, on top of which superconductivity 
emerges under doping, was proposed as a Resonating Va- 
lence Bond (RVB) state where localized electrons are 
paired up into singlet (bond) dimers due to antifcrro- 
magnetic (AF) exchange. Here, a substantial energy is 
gained from quantum resonance between different (short- 
range) dimer coverings. Re-formulating these ideas at 
the Hamiltonian level, Rokhsar and Kivelson constructed 
a Quantum Dimer Model (QDM) involving orthogonal 
dimer coverings and local dimer flips 1 (see Fig.nja)), 
aimed to capture the physics of systems that possess a 
spin (pseudo-) gap. Surprisingly, the QDM proved also 
to be relevant in a variety of distinct fields such as e.g. 
frustrated Ising models [ij , spin-orbital models [H or su- 
perconducting junction arrays @. Although the ground 
state (GS) of the QDM on the square lattice is a Valence 
Bond Crystal (VBC) breaking lattice symmetry jjjj], finite 
doping is expected to melt the crystal. In that respect, 
the doped QDM is of great interest since potentially more 
tractable than microscopic t— J or Hubbard models while 
effectively retaining their basic physical processes (spin 
exchange and hole hopping). 

The exact nature of holons, the charged spinless exci- 
tations of the doped RVB state, has been debated over 
the last two decades. In particular, their statistics is still 
unclear. Earlier work based on a variational RVB wave- 
function [|| suggested that hole excitations are (weakly 
interacting) fermions. Such a conclusion, although only 
justified for small enough kinetic energy, was reproduced 
in the context of the doped QDM [9(. Furthermore, it 
was argued Q that the holon statistics should in fact 
be dictated by energetics considerations and that, un- 
der some circonstancies, a holon could bind to a "vor- 
tex", (see e.g. Ref. d for a simple definition) leading 
to a bosonic composite. Zi gauge theories have recently 
brought powerful new tools [l0( to describe such phenom- 
ena. However, a more quantitative large-scale numerical 
investigation of the microscopic doped QDM is clearly 



needed. In this Letter, I perform such a program using 
Exact Diagonalizations (ED) of periodic 6x6 and 8x8 
clusters [ill ] and propose a semi-quantitative phase di- 
agram. In contrast to previous investigations of doped 
QDM's [IH with non-positive off-diagonal matrix ele- 
ments (named "Frobenius" QDM's), I consider here the 
non-Frobenius doped QDM which incorporates the Fermi 
statistics [l3| of the original electrons. Major differences 
in the phase diagrams shown in Fig. [2] are observed and 
discussed in the text. Of particular interest are: (a) a 
new unconventional d-wave superconducting phase (with 
possible lattice symmetry breaking) and (ii) a spectacu- 
lar change of the hole statistics, from Fermi to Bose, by 
increasing doping or kinetic energy. 
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FIG. 1: (color online) Dimer flip (a) and holon hopping (b) 
in the QDM. 



The non-Frobenius doped QDM model. — Following 
Ref. 0] and recent investigations 12 [, I start with the 
quantum hard-core dimer-gas on the two-dimensional 
square lattice defined by the Hamiltonian: 
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where the sum on (c) runs over all dimer coverings (con- 
taining a fixed amount of Nh vacant sites), N c is the 
number of fhppablc plaqucttes, the sum on (c',c) runs 
over all configurations |c) and |c') that differ by a sin- 
gle plaquette dimer flip, and the sum on (c",c) runs 
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over all configurations c) and \c") that differ by a sin- 
gle hole (or vacant site) hopping along a plaquette di- 
agonal as pictured in Fig. [1] In this formulation, bare 
holons (i.e the moving vacancies) have Bose statistics. 
Since the system is originally composed of electrons of 
Fermi statistics, a faithful description should then assume 
that dimer configurations (with vacancies) are created by 
sets of (spatially symmetric) dimer operators expressed 
in the fermionic representation as e.g. in Refs. @, 0], i.e. 

d h = Tf^iT^Ji + fftfii) where 4j creates an electron 
of spin a and i and j are neighboring sites. The origi- 
nal AF bond couplings yield then, in the QDM language, 
J < and hence a non-Frobenius Hamiltonian with pos- 
itive — J = | J | off-diagonal matrix elements. Note that 
a bosonic convention for the dimers can equally well be 
used as will be discussed later on. Sofar the sign of t is 
not specified [l5| . Phase diagrams (Fig. ^ are investi- 
gated as a function oiV/\J\ and \t/J\. 
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FIG. 2: (color online) Phase diagram of the doped QDM's 
versus doping (x) and V/\J\ (a,c) or \t/J\ (b,d). The Frobe- 
nius (a,b) and non-Frobenius (c,d) cases labelled on the plots 
as J > and J < 0, respectively, (see text) are compared. 
The x = insulator has VBC order. The symbols in (a,b) 
[resp. (d)] are obtained from the data of Fig. [3] [resp. Fig. [J. 



About Phase Separation — On finite clusters, Phase 
Separation (PS) is signaled by a negative compressibil- 
ity i.e. the local curvature of the energy vs doping e(x) 
curve. In that case, the equilibrium hole densities of the 
two-component mixture can be obtained via a standard 
Maxwell construction as for an ordinary liquid-gas first- 
order transition. Using such a procedure one shows that 
the Frobenius QDM (i.e. with J > 0) phase separates 
at small doping [!§]■ The ED results for both signs of 
J are compared in Fig. [3l while for J > the x — > 
curvature of e(x) is always negative, it is always positive 
for J < 0, even for small \t/J\ ratios. This shows that 
the two models behave quite differently in the immedi- 
ate vicinity of the x = axis as can already be noticed 
on their respective phase diagrams (Fig. ^ : (i) PS for 



J > (consistently with the data of Ref. [f2| obtained for 
smaller \t/J\ and closer to the RK point V/J = I) and 
(ii) a homogeneous phase (to be identified next) in the 
non-Frobenius case. In the latter case of most interest 
here, I find a very linear behavior of e(x) at finite doping 
evolving into a region of negative curvature as soon as 
\t/J\ > 0.5 suggesting that some form of PS might occur 
at finite doping. 
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FIG. 3: (color online) On-site GS energy (per site) vs doping 
calculated on a 6 x 6 cluster. Full (open) symbols are used 
for the non-Frobenius (Frobenius) QDM labelled as J < 
(J > 0). The linear behavior at x — > has been subtracted 
for convenience, (a) V = and different values of |i|; (b) 
\t/J\ — 0.5 and different values of V (all in units of |J| = 1). 
Insets: mean-slopes between x = and doping x plotted vs x. 
For J > 0, a polynomial fit gives a minimum at x c providing 
the range [0, x c ] of PS (Maxwell construction) reported in 
Fig. [2ja,b) . For J < no minimum at finite x is seen. 



Melting of VBC order— How does the x = VBC 
order evolve under finite doping is of great interest. Al- 
though, unfortunately, no finite-size scaling is possible 
here [13|, a comparison with the Frobenius model (for 
which large systems could be studied, see 12|) is very in- 
structive. As shown in Fig. IDJa), the relevant q = (tt, 0) 
dimcr-dimcr correlator, decreases rapidly with the hole 
kinetic energy. However, at small enough \t/ J\, the 
magnitude of the structure factor for the non-Frobenius 
model is always, on our clusters, larger than the same 
quantity in the Frobenius model. Since the later was 
shown to remain finite in the thermodynamic limit [l2T ]. 
this strongly suggests that there is, also in the non- 
Frobenius model, a finite region in the vicinity of x = 
where VBC order survives. Note however that this phase 
is not phase-separated as in the Frobenius case. Interest- 
ingly enough, Fig.[4|a) also reveals at larger \t/J\ a sud- 
den drop of the structure factor that we shall attribute 
later on to a real transition. 

The d-wave superconductor — Let us now refine the 
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characterization of the low doping phase in the non- 
Frobenius model. Two doped holes [3] are found to 
form a zero-momentum bound-state, of d x 2_ y 2 symme- 
try providing one chooses t > [15(. This is revealed 
by the "quasi-particle" (QP) peak at the bottom of the 
two- hole spectral functions plotted in Fig. [5{ a-f). The as- 
sociated weight Z2h vanishes when approaching the RK 
point, i.e. V/\ J\ —* 1, where VBC correlations become 
algebraic (Fig.|H[h)). This suggests that pairing is due to 
confinement by the VBC, at least at large length-scales. 
Incidently, a few higher energy peaks with approximate 
(tjt — l) 2 / 3 energy dependancc can be associated to so- 
called "string resonances" , typical of confinement phe- 
nomena. Interestingly enough, I also observed that the 
IVJl-dependance f Z 2h (Fig.^g)), Z 2h ~ |J/i|", v < 1, 
is quite similar to the case of a two-hole pair propagat- 
ing in a quantum AF [l6j |. Since PS docs not occur in 
the vicinity of the x = axis one then expects a d-wavc 
hole-pair superconductor in some finite doping range (see 
Fig. [U[c,d)). One could even argue (see above) that, at 
sufficiently small doping, d-wave (pair) superconducting 
and VBC orders are likely to coexist, although VBC or- 
der is not required for superconductivity. 
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FIG. 4: (color online) (a) Columnar (i.e. q = (tt,0)) VBC 
structure factor (normalized to its x = value) versus \t/J\. 
(b) "Fermionic" and "bosonic" signs (as defined in text) vs 
\t/J\. (a) and (b) correspond to the non-Frobenius doped 
QDM with V = (circles) and V/\t\ = 0.2 (lozenges) at 
x = 1/6 (computed on a 6 x 6 cluster). In (a) additional data 
for the Frobenius case are shown for comparison. Vertical 
dotted segments give the approximate separations between 
the phases discussed in the text. 



Holon statistics: Fermi vs Bose — Before turning to 
the investigation of the actual statistics of the holons, 
it is useful to mention briefly an alternative convention 
for the non-Frobenius QDM. Indeed, dimer operators 
can be represented in the usual bosonic representation 



4- -J 

operator b ia carrying spin a = hence leading to a 
change of the sign of J, i.e J > (provided is ori- 
ented e.g. from the A to the B sublattice). However, to 
preserve the fermionic character of the original electron, 
the bare holon should be a fermion, meaning that one 
now keeps track of some fixed (arbitrary) ordering of the 
holons [17| ■ Obviously, the actual statistics of the holon 
is independent on the representation used. 
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FIG. 5: (color online) (a-f) Spectral functions associated with 
the propagation of a q = (0, 0) d-wave pair of holes versus 
frequency ui (computed on a periodic 8x8 cluster) for fixed 
V = (left) or \t/J\ = 1 (right). The weight of the hole-pair 
QP peak (at the lowest energy edge) is plotted vs \J/t\ (g) 
and V/\ J\ (h). Comparison with data on a 6 x 6 cluster in (g) 
suggests weak finite size effects. Data for the t-J model (stars) 
from Ref. [TH are also shown for comparison in (g). Note that 
"string resonances" can be seen above the QP peak in (a-f). 



Motivated by [8( I now define "fermionic" (R=F) and 
"bosonic" (R=B) average signs as: 

si g n R = J2\J2 <**i c > i«i c >i i/E k*o ic>r , 

H a c£H a c 

where the first sum is performed over all classes H. a 
of configurations with fixed hole positions, the second 
sum runs over all dimer configurations within each class. 
The index "R" refers to one of the two (equivalent) rep- 
resentations of the non-Frobenius model (as discussed 
above) i.e. choosing bare holons with either Fermi (R=F) 
or Bose (R=B) statistics together with J > (R=F) 
or J < (R=B). Therefore, if holons truly behave as 
fermions (resp. holons) one expects Signp ~ 1 (resp. 
Sign B ~ 1). Since a fermion (resp. boson) can be seen as 
a bare bosonic (resp. fermionic) holon bound to a vortex 
(or fluxoid) [1, 0, H(| , one also expects simultaneously 
Sign B ~ (resp. Sign f ~ 0) as the presence of vortices 
will "mess up" the equal sign of the weights of configura- 
tions with fixed hole positions. From the results displayed 
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in Fig. 0|b) one clearly identifies a clear and simultane- 
ous rapid crossover of Sign F and Sign B signaling a change 
of statistics of the physical holons by binding/unbinding 
of vortices: at small (large) kinetic energy/doping, holon 
behaves as fermions (bosons) as argued from previous 
analysis [1, Q. The separation line between these two 
regions have been estimated and reported in the phase 
diagrams of Fig. [2] (c,d) as a thick line. At large kinetic 
energy one expects the holons to Bose condensate giv- 
ing rise to a charge-2e superfluid similar to the one of 
Refs. PEE! and Fig. [2] (a,b). The dotted lines of Fig. H 
(c,d) limit the area at small t or x where Sign F = 1 
within better than 5 x 10~ 3 , possibly connected to the 
supcrconduting and/or VBC phase. 

Discussion and conclusion — To summarize, a non- 
Frobenius doped QDM retaining the Fermi statistics 
of the original electrons reveals a rich phase diagram. 
Holons doped in the VBC host are shown to behave as 
fermions at small enough doping and to pair up. This 
leads to a d-wave superconductor (possibly coexisting 
with VBC order) which shows striking similarities with 
the 2D t-J model. A key feature that permits the ex- 
istence of this phase is the absence of PS in the imme- 
diate vicinity of the Mott insultor, in contrast to the 
Frobenius QDM [13]. At larger doping, with the disap- 
pearance of VBC order, holons behave as bosons and are 
expected to Bose condensate leading to a charge-2e su- 
perfluid [3, 12 1 . In this "statistics transmutation" the 
role of vortices (more specifically vortices as described 
in Ref. [h]]) is outlined. A complex behavior is also en- 
countered in-between characterized by a rapid cross-over 
region where Sign F ~ Sign B and some tendency towards 
PS. Identifying more precisely this cross-over is beyond 
the power of present-day computers but various hypoth- 
esis can be formulated. One scenario is a region of PS 
(i.e. macroscopic coexistence) of the two limiting phases 
which have been clearly identified in this work, as sug- 
gested by a slight negative curvature of e(x) vs x in Fig. [3] 
Another alternative is that both (charged) fermionic and 
bosonic excitations coexist at intermediate doping in a 
way which might bear some resemblance with the "holon- 
hole superconductor" of Ref. [HI with, w.r.t. the large- a; 
superfluid, additional low-energy excitations from paired 
holes. Incidentally, the observation of quantum oscil- 
lations attributed to some Fermi surface has been re- 
ported recently in cuprate superconductors [l9| . In any 
case, the tendency towards phase separation leaves room 
for modulated/unidirectional structures (or stripes) once 
long-range Coulomb repulsion is included. This is of cen- 
tral interest since modulated structures have been seen 
in cuprate superconductors, e.g. in recent STM experi- 
ments m. 

Note that the present study has been realized on the 
square lattice where the undoped QDM has VBC order. 
It would also be of great interest to investigate the case 
of the triangular lattice where a deconfincd RVB GS can 



be realized in the Mott insulator 21 1. 
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